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Abstract 

Interacting quantum systems evolving from an uncorrelated composite initial state generically 
develop quantum correlations — entanglement. As a consequence, a local description of interacting 
quantum system is impossible as a rule. A unitarily evolving (isolated) quantum system generi- 
cally develops extensive entanglement: the magnitude of the generated entanglement will increase 
without bounds with the effective Hilbert space dimension of the system. It is conceivable, that 
coupling of the interacting subsystems to local dephasing environments will restrict the genera- 
tion of entanglement to such extent, that the evolving composite system may be considered as 
approximately disentangled. This conjecture is addressed in the context of some common models 
of a bipartite system with linear and nonlinear interactions and local coupling to dephasing envi- 
ronments. Analytical and numerical results obtained imply that the conjecture is generally false. 
Open dynamics of the quantum correlations is compared to the corresponding evolution of the 
classical correlations and a qualitative difference is found. 

PACS numbers: 03.67.Mn,03.67.-a, 03.65.Ud, 03.65 Yz 
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I. INTRODUCTION 



The exploration of the nature and the extent of correlations generated by the many-body 
dynamics has both fundamental and practical applications. One of the fundamental issues 
in the investigation of many-body dynamics is finding an optimal set of coordinates [1, 2]. 
This problem is solved in classical mechanics by introducing partition of a complex system 
into smaller subsystems, i.e., introducing of degrees of freedom. Description of the composite 
system is furnished by local descriptions of the subsystems. The adequacy of a particular 
partition depends heavily on the nature and the extent of the correlations between the local 
degrees of freedom. 

The role played by correlations in classical and quantum mechanics is substantially differ- 
ent. This is due to the presence of the quantum correlations, or entanglement, in a composite 
quantum state, having no analog in the classical world [3]. In contrast to classical corre- 
lations [4], extensive entanglement makes the partition of a quantum system meaningless, 
since local measurements do not provide information on the state of an entangled system 
[3]. 

The problem of the optimal partition is deeply connected to the foundation of many-body 
dynamical simulations. The complete description of system composed of fully correlated sub- 
systems should grow exponentially with the number of subsystems involved. A possibility of 
representing a state of a complex system as a mixture of independently evolving uncorrelated 
states (trajectories) solves in principle the problem of many-body simulations, permitting to 
sample single trajectories for simulation and averaging the result subsequently [5, 6]. This 
possibility is inherent in classical mechanics but is nongeneric in quantum case, due to the 
fact that typical interaction of quantum system generates entanglement. If the growth of to- 
tal (i.e. quantum and classical) correlations becomes restricted, the quantum dynamics can 
be efficiently simulated [7-11]. Nonetheless, it is still an open question whether restrictions 
on quantum correlations alone are sufficient to provide for efficient simulations [12]. 

Addressing the problem of dynamical generation of correlations it is necessary to dis- 
tinguish between the unitary evolution of an isolated system and the open evolution of a 
system coupled to an environment. While a given unitary evolution can generate exten- 
sive entanglement, coupling the system to an environment is generally expected to restrict 
entanglement generation. This expectation originates in the general philosophy, seeing in 
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environmental-induced decoherence [13] the universal route of quantum-to-classical transi- 
tion. It is consistent with some established results on open-systems entanglement dynamics. 

Evolution of quantum correlations under the influence of environment was investigated 
both in the context of quantum to classical transition [13] and in the context of quantum 
information processing [14]. Most studies have been concerned with dynamics of entan- 
glement between noninteracting systems coupled to a bath. It was found that coupling 
to common environment is able to entangle noninteracting systems [15, 16]. On the other 
hand, coupling to certain local environments leads to total disentanglement of the systems 
in finite time [17-22]. The rates of disentanglement were calculated in bi- and multipartite 
systems of noninteracting qubits [19, 23-26] and quidits [27, 28], locally coupled to various 
environments. A number of studies addressed dynamical generation of correlations between 
interacting subsystems in the presence of the environment. Production of entanglement 
between qubits, modeling a system of ions, coupled to environment through their center 
of mass motion in ion-traps, was investigated in Ref. [28]. It was found that the coupling 
to environment diminishes the maximally achievable entanglement, with the corresponding 
entanglement loss increasing with the number of ions. Ref. [29] explored dynamics of entan- 
glement in quantum Heisenberg XY chain, immersed in a global purely dephasing bath. The 
robustness of entanglement against the dephasing was related to the number of spin in the 
chain. Coupling of interacting subsystems to local environments was considered in Refs. [30] 
and [18]. Ref. [30] investigated the generation and transfer of entanglement in harmonic 
chains. The creation of entanglement by suddenly switching on the interaction in the chain 
was found robust against the decoherence induced by coupling of the oscillators to local 
harmonic baths. The model of two harmonically coupled quantum Brownian particles was 
treated in Ref. [18]. It was found that in the physically interesting range of parameters the 
interaction between the particles cannot prevent their eventual disentanglement, induced by 
coupling to local baths. 

While the observed disentanglement of the noninteracting systems by coupling to local 
environments meets the common intuition about the quantum-to-classical transition, the 
picture of dynamics of entanglement in the presence of interaction is not so clear. Searching 
for an efficient and a universal environment-induced mechanism of restricting the extent 
of the generated entanglement, it seems necessary to focus on the following aspects of the 
dynamics of correlations. 
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First, the scaling of the generated correlations with the effective Hilbert space dimension 
of the interacting subsystems must be considered. This is in contrast to the context of 
the quantum information processing where the object of interest is usually the scaling of 
entanglement with the number of degrees of freedom (qubits). The expectation is that the 
environment-induced restriction on the generation of entanglement becomes most significant 
in the range of the large quantum numbers of the system, which is commonly associated 
with the quantum-to-classical transition. In fact, extensive entanglement in the large Hilbert 
space dimension seems impossible without creating the " cat-state" superpositions, which are 
expected to be destroyed by the decoherence. 

Second, the dynamics of correlations must be followed on the short, interaction time 
scales. It is possible that the long-time dynamics of an open composite system, approaching 
equilibrium, is disentangled, but the entanglement generated on the interaction time scales 
is so large that the partition of the system has no meaning. 

Moreover, since a common environment will generically entangle noninteracting systems, 
coupling to local environments seems necessary to provide for a generic route to a disentan- 
gled dynamics. 

The present study focuses on the investigation of environment-induced constraints on 
the dynamics of quantum and classical correlations in the open bipartite composite system. 
The system consists of two nonlinearly interacting harmonic oscillators, coupled to local 
purely dephasing baths. The distinction between the dephasing and pure dephasing has 
first appeared in the context of NMR [31]. Pure dephasing corresponds to loss of coherence 
in the energy representation. The two prototypes of underlying stochastic processes leading 
to dephasing are the Gaussian and the Poissonian processes [32]. R. Kubo based his line- 
shape theory [33] on the Gaussian model. The Kubo's model is the cornerstone of the 
condensed-matter spectroscopy. Recently exceptions to the Gaussian paradigm have been 
found experimentally [34-36] in ultrafast vibrational spectroscopy. The Poissonian model 
was shown to describe the dynamics adequately. Quantum Poissonian stochastic models have 
first appeared in the gas collision theory. They are also employed in the condense phase 
physics. For example, the Poissonian noise has been considered as a source of decoherence 
in quantum dots [37-39]. Due to the fundamental and the experimental relevance of the 
Gaussian and the Poissonian stochastic processes they were chosen as the source of the 
dephasing in the present study. 
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The models aim to explore dynamics of correlations in a composite system of coupled 
multilevel subsystems at large effective Hilbert space dimension. Examples of such systems 
include multimode molecular vibrations [40], linear and nonlinear quantum optics [41] and 
cold trapped atomic ions [42]. The primary goal is to locate a generic mechanism by which 
the decoherence keeps an interacting composite system "approximately disentangled" all 
along the evolution. Dynamics of quantum and classical correlations and their scaling with 
the effective Hilbert space dimension of the system are compared. 

The measures of quantum and total, i.e. quantum and classical, correlations are defined 
in Section II. Section III examines the issue of the generation of quantum correlations (en- 
tanglement) in the model problems. Section IV presents numerical results on dynamics of 
both quantum and classical correlations and Section V summarizes the conclusions. 

II. MEASURES OF CORRELATION 

The state of a bipartite system is uncorrelated if it can be described by the form 

Pab = Pa®Pb- (1) 

A general correlated state can be Schmidt-decomposed [43] (Cf. Appendix A ) in the 
Hilbert-Schmidt space leading to: 

N 

p ab = ]T c t k\ ® Bj, (2) 

i 

where the sets {A} and {B} of operators are orthonormal in the Hilbert-Schmidt spaces of 
systems a and b. 

The number of non vanishing coefficients Cj in the Schmidt decomposition of a vector 
in a abstract tensor-product Hilbert space is called the Schmidt rank of the vector. To 
avoid confusion in the following presentation the term HS-Schmidt rank (or just HS rank 
for brevity) is adopted for the Schmidt decomposition in the Hilbert-Schmidt (HS) space 
of operators, while retaining the term Schmidt rank for the Schmidt decomposition in the 
corresponding Hilbert (pure) state space. A HS rank is a natural measure of total correlations 
present in a mixed state p (Cf. Appendix A). 

A special subset of mixed states is the set of separable or classically correlated states [4] . 
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The state is separable if it can be cast into the following form 

N 

P»b = ^2PiPa® Pi ( 3 ) 

i 

where < pi < 1 and Y^i Pi = 1 and p a and p b are density operators defined on the Hilbert 
spaces of the subsystems a and b, respectively. Separable states are mixtures of uncorrelated 
states, which can be completely characterized by local measurements. Therefore, partition 
of a composite system into parts has a strong physical meaning. Such a partition is always 
possible for classical probability density distribution of a bipartite system [4, 44]. 

Pure correlated states are always entangled. The measure of pure state entanglement can 
be defined by its Schmidt rank [45]. Estimating the measure of mixed-state entanglement is 
a difficult conceptual and computational problem [45]. One can look for decomposition of p a b 
into a mixture of pure states that are least entangled on average. The average entanglement 
corresponding to such decomposition is a possible measure of the mixed state entanglements. 
Unfortunately, such measures are notoriously difficult to compute. 

An alternative computable measure of the bipartite mixed state entanglement is the 
negativity [46] defined as follows: 

\\p Ta \\ - 1 

W) = ^ , (4) 

where X = Tr"\/xt£ is the trace norm of an operator X and T a stands for the partial 
transposition with respect to the first subsystem. The partial transposition T a , with respect 
to subsystem a of a bipartite state p a b expanded in a local orthonormal basis as p a ] y = 
Y,Pij,ki \i) 01 ® \k) is defined as: 

The spectrum of the partially transposed density matrix is independent of the choice of local 
basis or on the choice of the subsystem with respect to which the partial transposition is 
performed. The negativity of the state equals the absolute value of the sum of the negative 
eigenvalues of the partially transposed state. By the Peres-Horodecki criterion [47, 48] the 
negativity vanishes in a separable state. On the other hand, vanishing of the negativity does 
not imply separability of the state in general [48]. 

Finite negativity is necessary and sufficient condition for the presence of entanglement in 
particular type of mixed states, the so called Schmidt-correlated states [49-51]. In this case 
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the negativity can be related to the structure of the density operator, which facilitates the 
evaluation of the entanglement. 

The Schmidt-correlated states have the following form 

P = ^Pmn \<t>m) (<t>n\ <8> \Xm) (Xn| , (6) 
mn 

where S x = {\(fi m )} m =i and H 2 = {\Xm)} m= i are local orthonormal bases. Eq.(6) implies 
that p = Y^iPi (V'il) where = J2 m c m 1 0m) <8> |Xm) 2 for ever Y h i- e - all pure states 
in the mixture share the same Schmidt bases (Cf. Appendix A ) Si and S 2 . It has been 
proved [52] that for Schmidt-correlated states 

m<n 

i.e., the negativity equals half the sum of absolute values of the off-diagonal elements of 
the density operator, written in a Si ® H 2 local tensor product basis. It follows that the 
negativity of entangled Schmidt-correlated states is finite [52]. 

The negativity can be related to the structure of the density operator. Consider the 
density operator (6) having the following quasi diagonal structure: 

P= ^2 Pmn ^ 

n\ IXm )<Xn|, (8) 

\m-n\<A 

with A <C k. The sum of the absolute values of the off-diagonal elements can be estimated 
as follows: 

n=m+A n=m+A 

^ ^ I Pmn | ^ ^ | Pmn | 1 ^ ^ ^ ^ I Pmn | 1 *^ ^ ^ ^ ^ \J PrnmPnn 

niy^n mn m n=m— A m n=m—A 

(9) 

n=m+A n=m+A n=m+A 

<E E «=±^ = iyj e A-.+5E E p»„<2A, 

m n=m— A m n=m—A m n=m~A 

where the first inequality follows from the positivity of the density operator and the second 
is the inequality of geometric and arithmetic means. Therefore, 

N{p) < A (10) 

in the state (8). Since the negativity of maximally entangled state (corresponding to \p mn \ = 
1/k in Eq.(6)) equals (k—l)/2, as follows from Eq.(7), the negativity of quasi diagonal density 



matrices is negligible compared to the maximally entangled state. It should be noted that 
the form (8) with A <C k of the density matrix does not constrain the magnitude of the 
classical correlations present in the state. For example, a strictly diagonal matrix p mn = 5 mn 
corresponds to a maximally (classically) correlated separable state. 

Schmidt correlated states appear naturally in a composite bipartite dynamics admitting 
particular conservations laws [52]. The models of open-system dynamics considered in the 
following sections belong to that class. As a consequence, the presence and extent of entan- 
glement in evolving composite systems can be related to the structure of the density matrix, 
which can be inferred on the basis of relatively general scaling considerations. 

III. DENSITY OPERATOR OF A BIPARTITE SYSTEM UNDER LOCAL PURE 
DEPHASING. 

A. General considerations 

The model of open system dynamics considered is described by: 

= (A + £ 2 )p + Jp, (11) 

where the generators of local nonunitary evolution are Cj = —i[B.j, •] — TjVj, j = 1,2 and 
X = — ry[H 12 , •] stands for the interaction superoperator. The operators Hj are local system 
Hamiltonians, the operator Hi 2 is the nonlocal (interaction) term in the composite system 
Hamiltonian and Vj denote local bath-dependent superoperators. Coupling constants i\ 2 
and 7 measure respectively the strength of coupling to the local environments and the 
strength of the interaction between the subsystems. 

As a reference, the open evolution of noninteracting subsystems ( 7 = 0) is considered 
first. In this case a local dephasing evolution of each separate system takes place: 

j t P = (C 1 + C 2 )p. (12) 

In many models of open evolution [13, 53-55] it is found that the evolving state undergoes 
decoherence, characterized by the decay of the off-diagonal elements of the density operator 
in particular basis of the robust states. Ideal robust states retain their purity notwithstanding 
the interaction with environment. Examples of such models include interaction with a purely 
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dephasing environment, singling out energy states as the robust basis, quantum Brownian 
motion and damped harmonic oscillator at zero temperature T = 0, which select the robust 
basis of coherent states. While the robust states basis is determined by the type of the bath 
and the system Hamiltonian, the time scales of the decoherence generally depend on the 
initial state as well. 

Let us assume that the local superoperators C\ and £ 2 in Eq.(12) single out local robust 
states bases Si and S 2 . A composite noninteracting system evolving according to Eq.(12) 
from an arbitrary initial state is expected to decohere in the tensor product basis: Si <g> S 2 . 
That means that an arbitrary initial state density matrix will eventually diagonalize in 
this basis. Switching on the interaction between subsystems causes a competition between 
entanglement generation and decoherence induced by the local baths. For sufficiently weak 
interaction viewing the evolving density operator in the unperturbed tensor product basis 
Si ®S 2 of local robust states is a good starting point. If the interaction perturbs only slightly 
the evolution of an off-diagonal matrix element, it will decay on an almost unperturbed 
decoherence time scale. 

To proceed with a more quantitative argument the concept of the effective Hilbert space 
7i e ff is helpful. Since the energy of the evolving system is finite, the evolution can be 
effectively restricted to a Hilbert space with finite dimension. This Hilbert space is termed 
the effective Hilbert space of the system. Let A be a spectral norm [56] of the interaction 
superoperator I restricted to the effective Hilbert-Schmidt space (i.e., the space of linear 
operators on 7i e //) and A be a spectral norm of the dissipator V = V\ + V 2 restricted to this 
space. A and A correspond to the shortest time scales of the evolution generated by the 1 
and T>, respectively. When A < A, the interaction timescale is slow compared to the shortest 
decoherence time scale. As a consequence, the evolution of certain matrix elements is only 
slightly perturbed by the interaction. In that case the perturbed dynamics of the matrix 
element will follow essentially the course of the decoherence. Therefore, a rough distinction 
can be made between the region of the density matrix dominated by the decoherence and 
the region dominated by the interaction. The border between the two regions is defined by 
the condition 

nj = 0{\-\ (13) 
where is the unperturbed decoherence time scale of a matrix element p^-, i, j G Si <g> S 2 . 
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In case that the decoherence-dominated regions of the density matrix are not populated 
initially, they will stay unpopulated in the course of the perturbed evolution. This property 
will shape the structure of the evolving density matrix. If the states are Schmidt-correlated 
states with local Schmidt-bases Si and S 2 , being the local robust states bases, the relation 
can be established between the structure of the matrix and the entanglement of the state 
as indicated in the previous section. Qualitatively, the larger is the decoherence-dominated 
region the smaller is the negativity of the state. 

The relative extent of the decoherence- and the interaction-dominated regions in a given 
dynamics generally depends on the initial state and, in particular , on the effective Hilbert- 
space dimension k of the system. As a consequence, different scenarios may be expected at 
asymptotically large k. The growing contribution of the interaction-dominated regions will 
generally imply extensive entanglement generation. On the other hand, if the relative size of 
the interaction-dominated regions becomes negligible at large k the entanglement generated 
by the open system dynamics may be negligible or even asymptotically independent on k. 
This possibility appeals to one who believes in the environmental- induced decoherence as a 
universal instrument of quantum to classical transition. An interesting question is the fate of 
the classical correlations in this scenario. While decoherence-dominated dynamics can turn 
extensively entangled initial state into extensively classically correlated state, it is not clear 
that decoherence-dominated dynamics can generate extensive classical correlations when 
quantum entanglement is negligible all along the evolution. Negligible total correlations 
seem nongeneric and do not correspond to the intuitive picture of a "really interacting" 
system. Therefore, a scenario of negligible quantum and extensive classical correlations 
matches best to a generic mechanism of quantum to classical transition. 

B. Model calculations 

The model calculations are used to illustrate and verify the general considerations 
presented above. The evolution of a bipartite system is studied according to Eq.(ll) 
§- t p = (£1 + C 2 )p + Ip, where (Cj = -i\Hj, •} - TjVj, j = 1, 2 and J = -ry[H 12 , •] 
) with two types of dissipators, corresponding to the Gaussian [57, 58] and the Poissonian 
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[35] purely dephasing models: 



H,, 



H,-,p (Gaussian) 

V ^ p= < L , - ( 14 ) 

e l< P H jpe l0H3 — p (Poissonian) 

These dissipators have the Lindblad form [59] of a generators of quantum dynamical semi- 
groups. The Gaussian and the Poissonian generators are the two examples explicitely men- 
sioned in the seminal paper by G. Lindblad [59]. 

The model Hamiltonian is a simplified version of a nonlinearly interacting multimode 
system. The local Hamiltonians EL,, j = 1,2 are chosen to be Hamiltonians of harmonic 
oscillators: H, = Uj&^&j, where a] and aj- are the creation and annihilation operators, respec- 
tively. Two general types of interaction are considered. The first (X4), termed band-limited 
interaction, is motivated by the stimulated Raman interaction between the translational 
modes of ions in cold traps [60-62] (for example, in Ref. [61] the effective interaction be- 
tween the modes is reduced to the band-limited operator exp[±i7r(aj.a y + a x a£)]). The sec- 
ond type of interaction (X B ) is motivated by weakly nonlinear interacting modes emerging 
in molecular vibrations [40] and in nonlinear optics. The typical example in the nonlinear 
optics is the second harmonic generation modeled by the interaction Hamiltonian of the 
form H = %(a 2 b"l" + aJ" 2 b) [41, 60]. In addition, the dynamics of the cold ion traps can be 
operated in the regime, where the effective interaction is well approximated by a weakly 
nonlinear coupling [60-62]. 

The two types of interaction are generated by: 



-i 7 [AlA 2 + A^Ax,.] (=1 A ) 
- ?7 [(4) s (a 2 ) r + (4) s (ai) r , •] s = 1, 2, ..; r = 1, 2, .., (= X B ) 



(15) 



where Aj is defined by its matrix elements in local energy basis: (Aj) mn = 5 m) n-i- The 
structure of Aj assures that I a is band-limited with the spectral norm A = 0( 7 ). 

The important property of the dynamics Eq.(ll), with local dephasing Eq.(14) and in- 
teraction Eq.(15) is conservation of a particular additive operator in each case. The first 
type of interaction, Z^, preserves the number operator N = ajai + a|a 2 : 2^(N) = 0, which 
is also preserved by the local generators: (C\ + /^(N) = 0. The second type of interaction, 
X B , preserves the generalized number operator N rs = rk\ki + saJ^: X B (N rs ) = 0, preserved 
by the local generators, as well: (C\ + C^i^rs) = 0. 
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Assume a pure uncorrelated initial state |-0(O)) = \k0) (written in the local energies 
basis). The state |?/>(0)) is an eigenstate of N with the eigenvalue k. As a consequence, the 
first type of the interaction, 2^, will drive the initial state into a mixture of eigenstates of N 
corresponding to the eigenvalue k: pit) = J2 mn c mn \m l m ) (n l n \ with l m = k — m. Thus, k 
determines the effective Hilbert space dimension of the system in this case: dim(7Y e //) = k. 
Since | -0(0) ) is also an eigenstate of N rs with the eigenvalue rk, the second type of interaction, 
X B , will take it into a mixture of eigenstates of N rs corresponding to the same eigenvalue 
rk: p(t) = J2 mn c mn \m l m ) (n l n \ with l m = ^{k — m). The number of initial excitations k 
of the first oscillator determines the effective Hilbert space dimension (in the strong sense) 
of the system in this case: dim(H e ff) = k/s. To summarize: 

~M I 7 \ / 7 I J l rn = k-m {1 A ) 

p(t) = c mn \m l m ) {nl n \< (16) 

[ l m = r -(k - m) (I B ) 

In both cases, the resulting mixed state is a Schmidt-correlated state with a time- 
independent Schmidt bases. This property permits evaluation of the negativity of p in 
each case from its structure, as indicated in Section II. 



1. Gaussian vs. Poisson pure dephasing bath 

The difference between the two types (14) of environments can be understood from com- 
paring the local evolutions of a single oscillator coupled to the bath of each type: 

^ p = — i[H, p] — TV p. (17) 

In the Gaussian case, the Eq.(17) in the energy representation becomes: 

Prim i^mnPnm ^^mnPnm) (^) 

with u rnn = cu m — u n , leading to the solution p nrn (t) = p nrn (<S)e~' lu)mnt ~ T ' u ' 2mnt . Thus the effect 
of the purely dephasing Gaussian bath is the "diagonalization" of the density matrix in 
the energy basis (the robust states basis for this model) on the time scale that varies for 
different matrix elements p nm and increases with the distance \m — n\ of the element from 
the diagonal. The shortest decoherence time scale corresponds to the largest distance from 
the diagonal and decreases with the growing effective Hilbert space dimension of the system. 



12 



°0 0.02 0.04 0.06 0.08 0.1 

time (in r~ 1 units) 

FIG. 1: (Color online) Purity, HS-participation number and negativity of the density operator of 
two harmonic oscillators, evolving under local purely dephasing Gaussian (solid lines) and Poisso- 
nian (dashed lines) environments (Eq.(14)). The initial state is a pure maximally correlated state 
\ip) = -^==£^ = q \n) \k — n) for k = 2,4, 12. The coupling parameter to the bath F = 1 in both 
cases. The frequencies of the oscillators u\ = 0J2 = 1- While the decay rates in the Gaussian case 
depend on initial state and increase with the effective Hilbert space dimension k, in the Poissonian 
case the rates are practically independent of k. 

In the Poissonian case, the Eq.(17) in the energy representation becomes: 

Pnm i^rnnPrim E(c ^ l)p nm , (19) 

leading to the solution p nm {t) = p„ m (0)e~^ mn ' +r ^ Wm ™ Apparently, the robust states 
basis is once again the energy basis, but the decoherence rates of the matrix elements are 
limited by |r(e~ iaJmn * - 1)| = 2F", independently of the initial state. 

This difference in properties of the Gaussian and the Poissonian environments will result 
in different dynamics of the correlations in the composite bipartite system dynamics Eq.(ll). 
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2. Local dephasing driven dynamics 



To gain insight on the effect of dephasing on the correlations the simplest bath-driven 
dynamics is studied first, in which the composite system Hamiltonian vanishes altogether, 
meaning that no entanglement is generated during the evolution. The corresponding Eq.(ll) 
transforms into 



d_ 
dt' 

d_ 
dt' 



p = — Tj ^Rj, Hj,p (Gaussian dephasing), (20) 



3=1,2 



: p = - ^ (e'^pe^ - pj (Poissonian dephasing). (21) 

3=1,2 

Since the dynamics preserves local energies the effective Hilbert space dimension of the 
evolving system is determined by the energy range of the initial state. The initial state of the 
form = ^== Ylt=o \n) \k — n) ( \n) is a local energy eigenstate) is a maximally entangled 
state. It corresponds to the effective Hilbert space spanned by the states {\n) \k — n)} k n=Q . 

The solution to Eq. (20), the Gaussian case, is found: 

p(t) = ——^e'^l^+^lrnJt \n) \k -n) (m\ (k - m\ (22) 

ran 

when the solution to Eq.(21), the Poissonian case, becomes: 

p{t) = T^—Y e -CTi[l-e-^-» *] + r 2 [l-e-*-i,-»*])t | n) ^ _ n) H {k _ ml (23) 

k + 1 ^-^ 



rnn 



Decoherence rates in the Gaussian case (22) are r m * = Yiuj\^ mn + r 2 w| jjrm < A g = 
{rio;^ + T 2 ul mn ) and generally increase without bounds with the effective 
Hilbert space dimension k. For example, taking H,- — u^-ajaj, T r = (ri^ + r 2 cj|)(m-n) 2 
is obtained, with maximal rate A g = t oi } = (TiLuf + r 2 cjf)A; 2 . In the Poissonian case (23) 
the decoherence rates are bounded: t~\ = Re {Tifl - e ~ iu ^ mn<t> } + r 2 [l - e - ^ 1 -™**]} < A p = 

2(ri + r 2 ). 

Note, that both solutions (22) and (23) are Schmidt-correlated states. Therefore, the 
corresponding negativities can be calculated from Eq. (7) : 



M(p(t)) = - 1 ^ y~] e ( riW Lnn+ r 2^, mn )t (Gaussian dephasing), (24) 

Af(p(t)) = -i- e-^ili-e-^—^I+rali-e- 4 " 1 —*])* ( PoiSSOman dephasing) . (25) 

m<n 
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From Eqs.(24) and (25) it follows that both types of the purely dephasing dynamics ( 
Eqs.(20) and (21)) lead eventually to a complete decay of the quantum correlations (note 
that since the evolving state is Schmidt-correlated, its negativity vanishes if and only if 
the state is disentangled [52]). But the dependence of the time scales of the decay on the 
effective Hilbert space dimension k is different in the two cases. In the Poissonian case the 
rate of the negativity (25) decay is bounded by A p = 2(Ti + T 2 ), independent of k, while in 
the Gaussian case (24), the bound is A g = max mi „ <fc {Tiw\ mn + r 2 w| mn }, which generally 
grows with k. 

The total correlations (and, as a consequence, the classical correlations) follow a differ- 
ent course of evolution. The HS-rank (and HS-participation number) of initial state is k 2 
(see Appendix A for calculation of HS-rank of a pure state). The stationary solution corre- 
sponding to both Eqs.(22) and (23) is p st = j^j J2 m \ m ) 1^ ~ m ) ( m \ ~~ m \ with HS-rank 
(and HS-participation number) equal to k. Therefore, although the total correlations decay 
in both models, the stationary solution contains extensive classical correlations, i.e. the 
correlations that grow without bounds with the effective Hilbert space dimension k. 

Fig. (1) displays the negativity, HS-participation number and purity of the composite 
state evolving under Gaussian (20) and Poissonian (21) dephasing dynamics, corresponding 
to IT = ^ajaj, Ti = T 2 and initial state of the form = ^== Ylt=o \ n ) \k — n). The 
effective Hilbert space dimension is varied: k = 4, 12. As anticipated from the difference 
of the two types of environments, the decay rates in the Gaussian case depend on the initial 
state and increases with the effective Hilbert space dimension, while in the Poissonian case 
the rates are effectively independent of the initial state. 

3. Full dynamics 

At this point the interaction between the oscillators are introduced and the full dynamics 
according to Eq.(ll) with 7 7^ is followed. We shall consider a pure uncorrelated initial 
state of the composite system: |-0(O)) = |/c0), i.e. the state corresponding to the excitation 
of the /c'th level of the first oscillator and the ground state of the second. In that 
shown above, each type of the interaction (15) and the dephasing (14) considered admits 
a particular additive conserved quantity (a generalized number operator), which defines 
the effective Hilbert space of the composite system for each k and is responsible for the 
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remarkable property of the evolving state: the density operator is a Schmidt-correlated 
state in a time-independent Schmidt bases: p(t) = Ylmn Pmnif) \ m ) ( n \ ® Km) (l n \ (Eq.(16)). 
The Schmidt bases S x = {\m}} and H 2 = {|^ m )} are the robust (local energies) bases of the 
corresponding local open systems (17), with the correspondence m <-> l m , determined by 
the particular conservation law, depending on the type of interaction. This property allows 
us to relate the structure of the evolving density operator to its negativity, as indicated 
in Section II. The relevant structure of the evolving density operator is determined by the 
relative size of the decoherence- and the interaction-dominated regions of the corresponding 
density matrix. This structure is investigated for each types of interaction and dephasing 
and for different effective Hilbert space dimensions of the system. 

The overview in the preceding section of the dynamics driven solely by the local dephasing 
reveals important difference between the two types of local environment with respect to the 
anticipated structure of the evolving density matrix. In the Poissonian case the decoherence 
rates are of the order of the system-bath coupling: r m l < 2(Ti + T 2 ), as shown above. 
Therefore, evolution of the matrix elements is dominated either by the decoherence or by the 
interaction depending on the relative strength of the coupling constants and independently of 
the effective Hilbert space dimension. In models with weak system-bath coupling, r 12 <C A, 
the structure of the evolving density operator will only slightly be effected by the coupling 
to the Poissonian bath on the interaction time-scale 7 _1 <C r mn . As a consequence, the 
quantum correlations will develop almost unperturbed on the interaction time scale. 

A different dynamical pattern is anticipated in the case of the Gaussian purely dephasing 
bath. The decoherence rates in this case are 



£ m ,„_i, and 1 B = -i7[(a|) s (a 2 ) r + (4) s (ai) r , •], with (%) m „ = y/m5 m , n -i, (see Eq.(15)). In 
each case, the decoherence rate increases with the "distance" \m — n\ from the diagonal. As 
a consequence, the evolving density operator obtains a quasi-diagonal structure in the local 
energies basis, with the width A of the interaction-dominated region about the diagonal 
depending on the type of interaction. 

Let us assume for simplicity that Tiuj\ = Y 2 oj\ = T. In that case a matrix element p mn 
decoheres on the time scale r mn = [2T(m — n) 2 ] -1 . The spectral norm of X A is A = 0(7). 




■j )mn 



(26) 
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Therefore, the width about the diagonal of the evolving density matrix can be estimated 
from Eq.(13) as 



a = o(v^7t), 



(27) 



where T 7 is assumed. The spectral norm of X# is A = Olk 2 ^), where k/s is the 
effective Hilbert space dimension of the system (see Eq.16). As a consequence, from Eq.(13) 
A becomes: 



(28) 



In the band limited interaction case (X4) the quasidiagonal structure of the density oper- 
ator emerges (Eq.(27)), while in the case of the nonlinear interaction (Xg) a quasidiagonal 
structure is expected only if the nonlinearilty is weak: s + r < 4 (Eq.(28)). 

Perturbation theory supports the scaling considerations. For a normalized eigenoperator 
6; of the local evolution generator C\ + C\: (C\ + £\)Oi = XiOi- The interaction X perturbs 
the evolution. The action of the perturbed generator on Oi gives: (4 + 4 + it)6, = 



A;(0« + If the trace norm of Si is small compared to unity: 



Si 



o, 



1, the 



evolution of 62 is only slightly perturbed on the time scale of (A;) -1 . Therefore, if Re[A^] < 0, 
the perturbed Oi will decay on a time scale of |Re[Aj]| to the leading order in Si . To each 
density matrix element p mn in the nonperturbed tensor-product basis of the local energy 
states (the robust states bases) there corresponds the normalized operator O mn = \ml m ) (nl n \ 
such that ^6 mr ^ = Tr |p6 mn | = p mn . Defining S mn by (C\+ C\+l r )6 rnn = A mn (6 m „+5 mn ) 
with \ mn = i(uji(m — n) + uo^ilm — l n )) — T[(m — n) 2 + (l m — l n ) 2 ] we obtain for the trace 
norm of S mn corresponding to the first type of interaction X4: 



o\l 



r y [f( m - n)] 2 + [f(l m - / n )] 2 + [(m - n) 2 + (l m - /„) 2 ] 2 
7 1 



< ° \T(m-ny + (l m -l n )\ 
and for the trace norm of S mn corresponding to the second type of interaction X#: 



(29) 



S„ 



o 



< o 



1 

T V [f(m - n)] 2 + If (l m - Z ft )] 2 + [(m - n) 2 + (Z, 
yjlr m m' + n% 



In) 



212 



T(m-ny + (l m -l n f 



(30) 
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The width A of the interaction-dominated region is estimated by solving 



7 1 

r (m - nf + (l m - Q 




(31) 



for the interaction generated by Ia and 



7 y 7 ^' + rggj = -. 

r(m- n y + (i m -i n y 

for the interaction generated by X B . Using Eq.(26), Eq.(31) is simplified to: 



(32) 




from which A 



2|m — n| = 



A/27/r is found, in compliance with the estimate (27), and 



Eq.(32) is simplified to: 



7 /r\ r / 2 \/[{k — m) r m s + n s (A; — n) r ] 



= 1. 



(34) 




In this case, the width about the diagonal A = 2\m — n\ will depend on m. The upper bound 
on A was calculated from Eq.(34) in two cases. First, for the linear coupling r = s = 1 gives 



results comply with the estimation Eq.(28). 

Figure (2) displays regions of the density matrix, dominated by the interaction, vs. 
regions, dominated by the decoherence, with the boundary between the regions deter- 
mined by Eqs.(33) and (34) for k = 10,20,40,50 and 7/T = 3. The figure repre- 
sents the composite system density matrices p = J2 mn c mn \m k — m) (n k — n\ and p = 
^ mn c mn \ m ^(k — m)} (n ^(k — n)\, corresponding to Eqs.(33) and (34), with m indexing 
the columns and n indexing the rows. The contours of Eqs. (34) are plotted for the linear 
coupling (r = s = 1) and the nonlinear coupling (r = 1, s = 2). The quasidiagonal struc- 
ture of the density operator is apparent. Both in the case of linear and nonlinear coupling 
between the oscillators, the width grows with the effective Hilbert space dimension. This 
is in contrast to the case of the band-limited interaction (Ia), where the width about the 
diagonal does not depend on the effective Hilbert space dimension k. 

To conclude, in contrast to the Poissonian type dephasing, in the Gaussian case the 
interaction-dominated regions are located about the diagonal of the density operator repre- 
sented in the local energies basis. Since the initial state \ip(0)) = \k0) corresponds to the 
density operator with an unpopulated decoherence-dominated region, this region will remain 




Second, for the nonlinear coupling r = 1, s = 2 gives A < y^/c 3 / 4 . Both 
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FIG. 2: (Color online) The density operator of the state, evolving according to Eqs. (11) for 
various interactions (15) and Gaussian type local baths (14) is represented in the product of lo- 
cal energies bases (Schmidt bases). Boundaries are indicated, in each case, separating the outer 
(off-diagonal) regions, dominated by the decoherence, from the inner (near diagonal) interaction- 
dominated regions. The interactions correspond to the band-limited case (case A, Eq.(15), dot- 
ted lines), linear coupling r = s = 1 (case B, Eq.(15), solid lines) and the nonlinear cou- 
pling r = 1, s = 2 (case B, Eq.(15), dashed lines). The density matrices in the band-limited 
case are of the form p = Y> mn c mn \ m k — m) {n k — n\ and in the linear and nonlinear cases 
p = T, mn c mn \m (r/s)(k — m)) (n (r/s)(k — n)\. The effective Hilbert space dimension corresponds 
to k = 10, 20, 30, 40, 50 in each case. See explanations in the text. 

unpopulated all along the evolution. As a consequence, the evolving density operator will 
stay in the quasidiagonal form. According to Eq.(10) the value of negativity is bounded by 
A in each case: J\f(p) < A. Asymptotically, i.e. as k 3> 1 for the band-limited interaction, 
\f~k ^> 1 for the linear interaction and \fk ^> 1 for the nonlinear case, the width about 
the diagonal becomes negligible compared to k. In this case, the generated entanglement 
is negligible compared to the maximal entanglement compatible with the effective Hilbert 
space dimension. 
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FIG. 3: (Color online)The negativity, the effective HS-rank and the HS-participation number of the 
evolving density operator: cases AG (solid lines) and A (dashed lines). In both cases u\ = u>2 = oj, 
T\ = T2 = I\ with Tuj 2 = (1/3)7 = (1/15)cj in case AG and T = in case A. Initial state 
\tp) = \k 0), with k = 4,6, ...,14. 

In the following section results of numerical calculations of the evolution of negativity, 
illustrating the foregoing discussion, are presented. The evolution of negativity is compared 
in each case with dynamics of the total (i.e., quantum and classical) correlations, as measured 
by the effective HS-rank and HS-participation number of the evolving density operator. 

IV. NUMERICAL RESULTS 

In the present section the results of numerical calculations of negativity A/"(p), HS- 
participation number k(p) and the effective HS-rank Xo.oi(p) (Cf. Appendix A) are dis- 
played and analyzed. The model is a bipartite composite state of two oscillators, evolving 
according to Eq.(ll). The dynamics simulated is classified according to the type of a local 
bath, Eq.(14), and the type of interaction, Eq.(15): 

• AG) The band-limited interaction Xa- Gaussian pure dephasing (Figs. (3,4)). 
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FIG. 4: (Color online) The negativity, the local energies and the HS-participation number of the 
density operator: cases AG and AP. Parameters: u)\ = u>2 = T± = T2 = T in both cases, 
Tuj 2 = 0.1257 = 0.025a; in AG (solid lines) and T = (1/15)7 = (l/75)w, <p = 2vr/7 in AP (dashed 
lines). Initial state \tp) = \k 0), with k = 4,6, 14. 

• AP) The band-limited interaction I a ■ Poissonian pure dephasing (Fig. (4)). 

• A) The band-limited interaction Xa ■ Isolated reference case (Fig. (3)). 

• BG1) The linear (r = s = 1) interaction X B . Gaussian pure dephasing (Figs. (5,6,7)). 

• BP1) The linear (r = s — 1) interaction Xb- Poissonian pure dephasing (Fig. (6)). 

• Bl) The linear (r = s = 1) interaction Xb- Isolated reference case (Figs. (5,7)). 

• BG2) The nonlinear (r = 1, s = 2) interaction Xb- Gaussian pure dephasing 
(Figs.(8,10,9)). 

• BP2) The nonlinear (r = l,s = 2) interaction Xb- Poissonian pure dephasing 
(Fig.(9)). 

• B2) The nonlinear (r = l,s = 2) interaction Xb- Isolated reference case (Figs. (8, 10)). 
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FIG. 5: (Color online) The negativity, the effective HS-rank and the HS-participation number 
of the evolving density operator: cases BG1 (solid lines) and Bl (dashed lines). In both cases: 
ui = u>2 = oj, in BG1 case: T± = T2 = T, Tu> 2 = (1/3)7 = (l/15)a> and in Bl case: ri = T2 = 0. 
Initial state = \k 0) for k = 4, 6, 24. 

In each case the evolution of the composite system starts from a pure uncorrelated state 
= \k0), where k is the initial number of excitations of the first oscillator, which determines 
the effective Hilbert space dimension of the system. 

Case AG (Figs. (3,4)). On Fig. (3) the negativity, HS-participation number and effective 
HS-rank of the evolving state in the presence of the bath is compared to the corresponding 
unitary evolution (case A). The amplitude of the negativity in the isolated case grows with- 
out bounds as the effective Hilbert space dimension k increases. Once the bath is introduced 
the amplitude of the negativity saturates to a value independent of k. On the other hand, 
both the HS-participation number and the effective HS-rank of the evolving state show that 
the total correlations grow without bounds when the effective Hilbert space dimension of 
the system increases. It is interesting to note the qualitative difference, most obvious in the 
unitary evolution (dashed lines), between the dynamics of the HS-participation number and 
the effective HS-rank on the shorter time scale, corresponding to the inverse frequency of the 
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FIG. 6: (Color online) The negativity, the local energies and the HS-participation number of the 
density operator: cases BG1 and BP1. Parameters: lo\ = uj2 = u>, T\ = T2 = T in both cases, 
Tlo 2 = (1/16)7 = (l/80)w in BG1 (solid lines) and T = (1/10)7 = (V 50 K <P = 27r / 7 in BP1 
(dashed lines). Initial state \tp) = \k 0), with k = 4,6, 18. 

oscillators cu -1 . While the HS-participation number is smooth on that scale, the effective 
HS-rank displays oscillations which follow closely after the corresponding dynamics of the 
negativity. 

Case AP Fig.4 compares the dynamics of correlations in case AP to case AG. The 
relative strength of couplings to different types of environments is chosen to match the time 
scales the local energies dephasing in both cases. It can be seen that in contrast to case AG 
both the negativity and the HS-participation number in case AP increase without bounds as 
the effective Hilbert space dimension grows similarly to the corresponding unitary evolution 
displayed on Fig. (3). 

Case BG1 (Fig. (5, 6, 7)). On Fig. (5) the negativity, HS-participation number and effective 
HS-rank of the evolving state in the presence of the bath is compared to the corresponding 
unitary evolution (case Bl). The amplitude of the negativity in the isolated case grows 
without bounds as the effective Hilbert space dimension k increases. In the bath-on case the 
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FIG. 7: (Color online) Negativity and HS-participation number (raised to various powers to fit 
linear dependence) measured at the first (in time) maximum on Fig. (5) for BG1 (solid line) and 
Bl (dashed lines) cases and on Fig. (6) for BP1 (dotted lines) case as a function of the number 
of excitations k (equal to dim(7i e /f)). Negativity in the BG1 case is raised to power 4 and the 
negativity in both the Bl and BP1 cases is raised to power 2. The powers of the HS-participation 
number are 3/2 in BG1 case and 1 in Bl and BP1 cases. 

amplitude of the negativity is obviously restricted but the quantita Both HS-participation 
number and effective HS-rank display the growth of the total correlations without bounds 
with the effective Hilbert space dimension of the system. Note the qualitative difference in 
dynamics of the two measures. 

Fig. (7) displays the maximal values of the negativity and the HS-participation number 
obtained in cases BG1, Bl and BP1 as functions of the effective Hilbert space dimension 
k. It is seen that in the Bl case the squared negativity and the HS-participation number 
scale linearly with k in compliance with the calculation in Appendix B. The same scaling is 
found in the case BP1. On the other hand, the negativity in the BG1 case scales as a fourth 
root of the effective Hilbert space dimension. The corresponding HS-participation number 
measuring the total correlations scales as k 2 ^ 3 . 

Case BP1 Fig. (6) compares the dynamics of correlations in case BP1 to case BG1. The 
relative strength of couplings to different types of environments is chosen to match the local 
energies dephasing rates. From this figure and Fig. (7) it can be seen that in contrast to case 
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FIG. 8: (Color online)The negativity, the effective HS-rank and the HS-participation number 
of the evolving density operator: case BG2 (solid lines) and B2 (dashed lines). In both cases: 
2uj\ = u>2 = ui, in case BG2: T\uj\ = I^c^l = (1/3)7 = (l/30)o> and in case B2: T\ = T2 = 0. 
Initial state = \k 0) for k = 4, 6, 28. 

BG1 both the negativity and the HS-participation number in case BP1 follow a dynamical 
pattern identical to the corresponding unitary evolution. 

Case BG2 (Fig. (8, 9, 10)). In Fig. (8) the negativity, HS-participation number and effec- 
tive HS-rank of the evolving state in the presence of the bath is compared to the correspond- 
ing unitary evolution (case B2). For a nonlinear interaction both the amplitudes and time 
scales of the dynamics depend on the initial state. As a consequence, the pattern of behavior 
changes with the effective Hilbert space dimension. This makes it difficult to compare the 
evolutions corresponding to different k. Comparing the open to the closed unitary evolutions 
for a fixed k it is seen that the global dynamics of the negativity is much stronger effected 
by the bath than the dynamics of the total correlations. For example, the total correlations 
may grow in the open dynamics similarly to the unitary case, while the negativity at this 
time is decaying in sharp contrast to the corresponding unitary behavior. 
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FIG. 9: (Color online) The negativity, the local energies and the HS-participation number of the 
density operator: cases BG2 (solid lines) and BP2 (dashed lines). Parameters: 2u>i = <jJi = uo in 
both cases, Yiu\ = r 2 w| = (1/8)7 = (l/80)w in BG2 and Ti = T 2 = (1/4)7 = (l/40)o;, <t> = 2vr/7 
in BP2 . Initial state = \k 0), with k = 4,6, ...,20. 

A possible way to compare values of the negativity and the total correlations at different k 
is to measure the values observed at the first maximum in the evolutions of these quantities. 
These measurement are displayed on Fig. (10). To understand the scaling, the negativity 
(squared to fit the linear dependence) and the HS-participation number obtained in cases 
BG2, B2 and BP2 are plotted as a functions of the effective Hilbert space dimension k. It 
is found that the negativity scales with \fk while the HS-participation number scale linearly 
with k. 

Case BP2 Fig. (9) compares the dynamics of correlations in case BP2 to case BG2. The 
relative strength of couplings to different types of environments is chosen to match the local 
energies dephasing rates. The negativity and the HS-participation number in case BP2 
follow a dynamical pattern identical to the corresponding unitary evolution displayed on 
Fig. (8). See also Fig. (10). 
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FIG. 10: (Color online) Negativity (squared to fit linear dependence) and HS-participation number 
measured at the first (in time) maximum on Fig. (8) for BG2 (solid line) and B2 (dashed lines) 
cases and on Fig. (9) for BP2 (dotted lines) case as a function of the number of excitations k (equal 
to dim(H e //)). 

V. SUMMARY AND CONCLUSIONS. 

Variety of open interacting bipartite systems were investigated in order to characterize 
restrictions, imposed by coupling to local environments, on the generation of classical and 
quantum correlations. 

The extent of the generated quantum correlations is determined by the interplay of two 
competing forces: the interaction, leading to development of entanglement, and the local 
decoherence, inducing a decay of entanglement. The relative magnitudes of the local deco- 
herence rates and the cut-off frequency of the interaction in the effective Hilbert space of 
the composite system determines the relative size of decoherence- and interaction-dominated 
regions of the density operator in local robust states basis. The presence of the decoherence- 
dominated regions constrains the structure of the evolving composite density operator, re- 
stricting the extent of entanglement, generated by the interaction. 

The character of restriction depends on the type of bath and the type of the interaction. 
The two different praradigms of the dephasing, the Poissonian and the Gaussian, lead to 
very different correlation dynamics. In models with band-limited decoherence such as the 
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Poissonian pure dephasing model, either the decoherence or the interaction dominates the 
dynamics, depending on the relataive strength of the coupling constants and irrespectively 
of initial state. Numerical calculations performed on a bipartite system of two interacting 
harmonic oscillators, coupled to local Poissonian baths, support this conclusion. 

Open systems with Gaussian pure dephasing belong to a different class of models. This 
class is characterized by unbounded growth of the decoherence time scales with the effective 
Hilbert space dimension of the system. As a consequence, constrains on the structure of 
evolving state and restriction on the extent of entanglement are generally expected. Still the 
precise character of the restriction depends on the type of interaction between the subsys- 
tems. Coupling local Gaussian environments to subsystems with band-limited interaction 
between them, imposes an upper bound on the extent of generated entanglement, which 
is independent of the effective Hilbert space dimension of the system. As a consequence, 
asymptotically, i.e. at sufficiently large effective dimension, the generated entanglement 
is negligible, compared to entanglement generated in the corresponding unitary dynam- 
ics. Interactions which are not band-limited generally produce extensive entanglement, 
notwithstanding the type of local environment. Nonetheless, in models with local Gaus- 
sian environments the scaling of entanglement with the effective dimension is limited by 
the local decoherence. The precise limit depends on the nonlinearity of the interaction. In 
the model of two nonlinearly interacting harmonic oscillators stronger nonlinearity implies 
weaker bounds on the generated entanglement. When the nonlinearity exceeds some maxi- 
mal value no restriction on the extent of entanglement is expected. Numerical calculations 
support these predictions. 

Estimation of bounds on negativity in the evolving state was based on analysis of the 
structure of the density matrix in particular local robust states bases. Relating the nega- 
tivity to the structure of the density operator was facilitated by the observation that the 
evolving states are Schmidt-correlated due to particular conservation laws observed by the 
interactions. The corresponding Schmidt bases are built of local robust states selected by 
local purely dephasing environments — the local energy bases. Since the presence of ex- 
act conservation laws is nongeneric in physical models, it should be noted that numerical 
evidence shows that the qualitative picture presented above is robust. 

Dynamics of the total correlations was investigated numerically to compare with the 
corresponding dynamics of the entanglement. It was found that evolution of the total (and, 
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as a consequence, classical) correlations display a different dynamical pattern. In the band- 
limited interaction model, the amplitude of the total correlations grows without bounds 
with the effective Hilbert space dimension, while the negativity tends to an asymptotic 
behavior independent of the effective dimension. In the linear interaction model, though 
the amplitudes of both the quantum and the total correlations grow without bounds with 
the effective Hilbert space dimension, the total correlations scale with a higher power of 
the dimension. In the nonlinear interaction, comparison is impeded by the fact that the 
evolution of both the entanglement and the total correlations display a variety of time- 
scales. Nonetheless, inspection of the numerical evidence shows that the total correlations 
always scale with a higher power of the effective Hilbert space dimension. These findings 
can be informally interpreted as a trade-off between the classical and quantum correlations: 
since the total correlations are (relatively) unaffected by the environment, restriction on the 
entanglement generation must be "compensated" by the growth of the classical correlations. 

Considering the restriction on the generation of entanglement, a natural question arises: 
is the observed restriction substantial, i.e. is the given partition of the composite system 
meaningful? When can a composite systems be regarded as approximately disentangled"! 
The answer depends on the definition of the relevant scale of a measure of entanglement in 
the evolving system. Is the scale unity or some power of the effective Hilbert space dimension 
or neither? 

One possibility is to compare the entanglement, generated in the open evolution to the 
entanglement, generated in the corresponding unitary evolution. Numerical evidence ob- 
tained in the present study shows that entanglement is always relatively restricted in the 
open system dynamics. In some cases, such as the Gaussian pure dephasing, it can even be 
negligible in asymptotically large Hilbert space dimensions. This comparison elucidates the 
role of the decoherence in constraining the generation of the quantum correlations. Never- 
theless, the magnitude of the entanglement generated in a particular open evolution may 
still be large in some absolute sense. 

An alternative scale of entanglement is set by the maximal entanglement compatible 
with the effective Hilbert space dimension. The results of the present study show that in 
some models, such as the Gaussian pure dephasing and weakly nonlinear or band-limited 
interactions, coupling to local environments does the job, i.e. it restricts the generated 
entanglement to bounds, negligible compared to the maximal compatible entanglement. 
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Still, in all cases apart from a band-limited type of the interaction, entanglement, generated 
on the interaction time scale in the open system evolution, grows without bounds with the 
effective Hilbert space dimension. As a consequence, this scale may become irrelevant in 
large effective Hilbert dimensions, due to a limited experimental resolution. 

To conclude, common models of local decoherence do not provide a universal pathway to 
an approximately disentangled evolution of a bipartite composite system in the presence of 
interaction. It follows that, contrary to expectations, coupling to local environments does 
not generally validate partition of composite quantum systems. 

APPENDIX A: THE SCHMIDT RANK AND THE HS-SCHMIDT RANK. 

The definition of the Schmidt rank of the bipartite composite state [3, 43] is reviewed. 
Let be a state in the composite Hilbert space 1~Lyi = Ji\ <8> T~ti- There exist a following 
representation (a Schmidt decomposition ) of the state: \ip) = J2i c i lOi ® 1^2' where \i) 12 
is an ortho normal basis in the Hilbert space 7^1,2- While a Schmidt decomposition is not 
unique, the set of non vanishing coefficients q is invariant (modulo irrelevant phases) under 
the local unitary transformations and is characteristic of the state This set is shown 
to be the square root of the spectrum of the reduced density operator of either subsystem. 
The number of non vanishing coefficients q is called the Schmidt rank xi^) °f the state and 
equals the rank of the reduced density operator of either subsystem: xijP) — ran k{Tri{pi 2 }}. 
To calculate the Schmidt rank of a state expressed in an arbitrary tensor product basis 
1 1.7)2 one calculates the rank of the matrix p 2 — C^C, where — c^-: 

Tr!{pi 2 } = E{rw,M} C V C lAnSkn \j) {l\ = £{iJ,J} %4 \j) ( l \ = E{j,o(^)/j \j) 

The Schmidt rank characterizes the extent of correlations present in the state. The 
uncorrelated (product) state has x = \ but generally xijP) < min {dim(7ii), dim(H2)}. 
The maximally correlated state has = min {dim(7ii), dim(7Y2)} and q = Cj, 

Generally, some of the coefficients q are much smaller than others and as a consequence 
dropping the corresponding contributions to the Schmidt decomposition does not lead to an 
observable effect. This suggests a definition of the physically reasonable effective Schmidt 
rank [9] Xe : XeijP) = xijP')) with = J2iei e c « N)i ® |*) 2 j where I e is the smallest set of 
indices such that |||^) — < e. An alternative measure is a participation number [63] 
= l/Trlpl 1 } with p 2 = Tri{pi 2 }. The participation number of a state, characterized 
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by M equal substantial contributions to its Schmidt decomposition is seen to be M, which 
motivates the definition. 

A mixed state displays both quantum (entanglement) and classical correlations. The 
extent of the total correlations can be characterized by the Schmidt rank of a density op- 
erator. With a slight abuse of terminology the term HS-Schmidt rank (HS indicating the 
Hilbert-Schmidt space) or just HS-rank is adopted. The definition of the HS-rank views 
the density operator of a composite system as a (unnormalized) pure state ("superket" [10]) 
in the Hilbert-Schmidt space of system operators. The Schmidt rank of the corresponding 
"superket" defines the HS-rank (denoted x(p)) °f the density operator. The notions of the 
effective Schmidt rank Xe(V') and the participation number can be transfered to the 
HS-rank of the density operator. For brevity, the corresponding measures of the total cor- 
relations are termed effective HS-rank and HS-participation number and denoted by Xe{p) 
and k(p), respectively. 

The calculation of the HS-rank proceeds as follows. Let p\ 2 = j k i} Pijki\ij) 
be a density operator of the composite system. In the superket notation it has the form 
\p) 12 = Y.{i,j,k,i} Pijki\\ij) The corresponding density superoperator is TZ 12 (p) = 

J2{i,j,k,i,i',j',k',i'} PijkiPi'fk'v llu) (^1) (\i'f) (k'l'W and the reduced density superoperator is 

K 2 {p) = Tn{^i 2 (p)} = ,j,k,l,i',f,k',l',m,n} PijklP*i'rk'l' S rniSnkS m i'Snk'\\j) (l\) (\f) (l'\\ = 

^{j,i,f, i'} R jW lb') ( l D (\f) ( l % where R jii'V = ^ikPijkiPi'j'k'V The HS-rank x(p) of the 
density operator p is x{p) = rank{7£ 2 (p)}. The effective HS-rank and the HS-participation 
number are calculated similarly. 

Finally, note that x(\i/>) = xC^) 2 - 111 fact, X% = (rank{p 2 }) 2 = rank{p 2 <S> 

pl) = T ^n2 { ij,i,kj',i'} a ij a u a kj' a ki' \j) ( l \ ® 10 U'\} = rank {T / {i, j ,i,kj',i'}pijki'p*ikj' \j) (A ® 

|0 = rank{£ {j tJ ,, v} Rjvif \j) (A ® 10 0*1} = rank{£ w i?^ \j) (f\ ® 10 (I'D = 
rank^^R^, \jl) (fl'\} = x(|^> M). 

APPENDIX B: CALCULATION OF THE EFFECTIVE SCHMIDT RANK OF 
THE COMPOSITE STATE OF TWO LINEARLY INTERACTING HARMONIC 
OSCILLATORS 

A system of two linearly interacting harmonic oscillators is considered with the Hamil- 
tonian H = cj(a|ax + a) 2 d, 2 ) + 7(a|a 2 + aJ,ai) The initial state is |-0(O)) = \0k) in the 
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local energies basis. The state at t > becomes: \ip(t)) = Ylt=o ° n 1^ — n n )' wnere 
c n (t) = y / fc!c ° s(7 ti( fc l n ff - 1! ^e-^. The width A fe of the distribution of expansion coef- 
ficients \c n \ 2 is estimated at t = n/A^y for fc 3> 1. This width is a reasonable estimate for the 
amplitude of the effective Schmidt rank of the state: ~ 

The distribution of the coefficients \c n (ir / A'j)] 2 = 2 k n \(k-ny. * s P ea ked around n = fc/2. 
To estimate A it is assumed that A fc <C fc. A fc is defined by: ^ | c n (71-/47) | 2 | n=n * = 0, 
where n* = fc/2 — A&/2. Performing the derivation under the Stirling approximation for 
the factorials (valid at fc 3> 1) leads to n *^- n *) — m2 (^?~)- For highly peaked distribution 

« 1, therefore ln 2 (^) « (^l 11 ) 2 . Also . * « | to the leading order in 
Finally f » (^) 2 » (^) 2 = (f)|) 2 from which A, = and X (^) « A fc = y/k. As 
follows from the relation x(\ip) (ip\) = xijPYi proved in Appendix A, the amplitude of the 
effective HS-rank scales as fc. 

The obtained result can be used to estimate the amplitude of the negativity in the pure 
state evolution. In fact Af(\^(t)) (^(t)\) = |(| E„ c n| 2 - 1) by Ref.[46]. Taking c n = ^= = 
tt= for the purpose of scaling we obtain M(\iP(tx / A^/)) (■0(7r/4'y) |) = \\fk~ for the amplitude 
of the negativity. 
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